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Abstract :- Aerodynamic noise prediction has been an important and challenging research area 
since James Lighthill first introduced his Acoustic Analogy Approach over fifty years ago. This 
talk attempts to provide a unified framework for the subsequent theoretical developments in this 
field. It assumes that there is no single approach that is optimal in all situations and uses the 
framework as a basis for discussing the strengths weaknesses of the various approaches to this 
topic. But the emphasis here will be on the important problem of predicting the noise from high 
speed air jets. Specific results will presented for round jets in the 0.5 to 1.4 Mach number range 
and compared with experimental data taken on the Glenn SHAR rig. It is demonstrated that non- 
parallel mean flow effects play an important role in predicting the noise at the supersonic Mach 
numbers. The results explain the failure of previous attempts based on the parallel flow Lilley 
model (which has served as the foundation for most jet noise analyses during past two decades). 
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1 Introduction 

Aeroacoustics is concerned with sound 
generated by forces and stresses resulting 
from the motion of a fluid rather than the 
externally applied forces and motions of 
classical acoustics. Music (i.e., harmonious 
sound) is often generated in this fashion. But 
this paper is concerned with aerodynamic 
noise, (i.e., cacophonous sound), which is of 
engineering interest because of the demand 
for quieter transportation vehicles-such as 
aircraft, busses and automobiles. In 
particular the noise from jet engine exhausts 
has been a major concern since the type of 
engine was introduced over fifty years ago 
and wind noise from automobiles is now 
receiving considerable attention since other 
previously dominant noise sources have 
been reduced to relatively low levels. 

The focus of this paper is on the 
fundamental theory. But since aerodynamic 
sound is just a byproduct of fluid motion 


and since the governing equations of that 
motion (i.e., the Navier-Stokes equations) 
are well established and have been for some 
time now (over 150 years) you might think 
that the only thing that we need to do here is 
write down these equations and be done 
with it. But the solution to these equations — 
especially for the turbulent flows that are the 
usual source of aerodynamic noise- has, to 
put it mildly, proved to be very elusive. In 
fact it is probably no exaggeration to say 
that fluid turbulence is still ranked among 
the great unsolved problems of classical 
physics. 

So there is a need to obtain useful and 
interesting information about the sound field 
without actually solving the full Navier- 
Stokes equations — which typically involves 
introducing some sort of reduced order 
model for these equations — which, in the 
present context, is usually referred to as an 
acoustic analogy (following James Lighthill 
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who was the first to do such analysis)[l]. 
This is, of course, analogous to the 
situation in turbulence modeling, which 
is almost always based on some fonn of 
the filtered (usually Favre [2] filtered) 
Navier-Stokes equations such as the 
Reynolds averaged Navier-Stokes (RANS) 
equations [3]. While there is some 
disagreement about the proper choice of the 
turbulence modeling equations, the situation 
in Aeroacoustics is considerably more 
contentious with significant disagreement 
about the appropriate starting equations 
(perhaps because there is no single set of 
equations that is optimal in all 
situations). There does, however, seem to be 
a consensus about some of the requirements 
for such equations. First of all, they should 
be derivable from the Navier-Stokes 
equations and secondly; they should, as 
argued in the well known text by Ann 
Dowling & Ffowcs Williams, be formally 
linear [4]. 


computational method, which means that it 
can, at best, be expected to provide an 
adequate representation of the most 
energetic part of the flow and perhaps (with 
enough computational effort) even the low 
frequency components of the sound field-but 
certainly not all of the spectral components 
that are of interest the high Reynolds 
number flows of practical importance. For 
example the noise from full jet engine 
exhausts has a peak frequency in the 200- 
500 Flz range, but the ear is most sensitive to 
sound in the 2,000 to 5,000 Flz frequency 
range. So jet noise prediction methods need 
to maintain their accuracy over an 
enormously broad range of frequencies. 

So as a practical matter the base flow 
will have to be determined from some 
reduced order model of the Navier-Stokes 
equation. To my knowledge, these models 
are almost always based on the usual 
hyperbolic conservation laws that can be 
written (fairly compactly) as [5, 6, 7, 8] 


2 Fundamental Equations 

It, therefore, makes sense to begin by 
dividing the flow variables in the Navier- 
Stokes equations, the density p, the pressure 

p, the enthalpy /?, the velocity v ( . , etc, into, 
say, their base flow components , 
p,p,h, and v. and ‘residual’ components 
p' , //, h' and , v \ , which are essentially 
defined by 

p = p + p', p = p + p', 

(1) 

h = h + h', v i =v i +v' i , 

with the idea being that base flow is to be 
found from some relatively inexpensive 


D oP = °> 


( 2 ) 


dp d ~ 
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(4) 


where the Latin indices range from 1 to 3, 
the summation convention is being used. 
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(7) 


D °f J i + ^M 


(5) 


for any function f , we assume that the base 
flow variables, as well as the original Navier 
-Stokes variables satisfy an ideal gas law 
equation of state, with y being the specific 

heat ratio, p e denotes a pressure-like 
variable that can differ from the 
thermodynamic pressure, and 0 . . with 

A = 1,2, 3, 4 denotes an, as yet, arbitrary 
4x3 dimensional stress tensor. (Greek 
indices will always range from 1 to 4.) 

These equations include, among other 
things, the Euler equations, the Navier 
Stokes equations themselves (which 
correspond to setting the first three 

components of 0 equal to the Newtonian 

stress and the fourth component to the heat 
flux vector). And more generally, the 4x3 
dimensional stress tensor can be inputted 
through a turbulence model, which also 
determines the difference between the 
effective pressure p e and the 
thermodynamic pressure p —but in a purely 
passive fashion. However, this difference is 
largely irrelevant since these equations will 
form a closed system of 5 equations in 5 
unknowns which can be solved for p e 

independently of its relation to the 
thermodynamic pressure p once a particular 
turbulence model has been introduced. 

It can now be shown that the remaining 
residual variables are determined by the five 
formally linear equations [5, 6, 7, 8] 


„ ov i 8 , 

D o u l +u J y— + y^Pe 
OX; OX; 


p' d ~ d „ 
r. = e 

p dxj " 8xj v 


8 ~ 

D ° P ' e + ~dx^ c2Uj + 
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r p ,dvj u ,df^ 


V 


dx ■ p dx 
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d " , f „ dv t 


( 8 ) 


where 


c 2 = y Rf -y^r 
P 


(9) 


is the base flow sound speed squared and 

T/ ^ S lJ Pe -^/ ( 10 ) 


denotes the total base flow stress tensor. 

These equations have been put into the 
linearized conservation law form by 
introducing the new dependent variables 

Pe=P-Pe+^ pv ' 2 (11) 


u i =P V i» 


( 12 ) 


as well as the new source strengths 
e”j , v = 1,2, 3, 4 given by the difference 


e m ~ e jv e jv 


(13) 


D °p' + -^- Uj =0, 


( 6 ) 


between the generalized Reynolds stress 
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e Y/ = 


-P V v Vy + 


Y- 1 ^ — 


8 V/ Pv'- +a v/ 


(14) 


y? 


(based on the residual velocities and 
enthalpy) and the base flow stress 6 • , 

where 


K = (r-i)K =( y - ! ) 


Y 1 A 

h' + -v' 2 
{ 2 ) 


(15) 


a 4i = — (Y-l)(^f ~°ij V 'j ) ( 16 ) 

with a ; y being the viscous stress, qi the 

heat flux vector and 8 vj denoting the 4x3 

dimensional Kronecker delta in the usual 
notation. 

The true non-linearity of these equations, 
which can be written more compactly in 
operator form, as 


= D viKj for 
p, v = 1,2,3, 4,5 


(17) 


where L^ IV and /3‘T arc first order linear 

operators and u v denotes the five 
dimensional solution vector 

{ M v} = { u i>P' e >P'}> (18) 

is hidden in the non-linear dependent 
variables (1 l)and (12) as well as in the non- 
linear source strength (14). The former 
non-linearity is again largely irrelevant here 
because the variable of principle interest (the 
pressure-like variable, p ' e ) reduces to the 

ordinary pressure fluctuation in the far field 
where the sound field is to be calculated. 

And in the acoustic analogy approach the 


latter non-linearity is dealt with by intruding 
a specific model for the generalized stress 
tensor (14). This is analogous to the 
situation in conventional turbulence 
modeling. But the modeling is now 
considerably more complex and being able 
to simplify the requirements by lumping all 
the modeling into a single tensor is an 
important consideration. So we might thi nk 
of requiring that the fundamental acoustic 
equations satisfy the following two 
conditions in addition to the two mentioned 
in the introduction: 

(1) The base flow about which the implied 
linrearization is carried out satisfies the 
usual hyperbolic conservation laws. And 

(2) A11 of the modeling requirements can be 

inputted through a single stress tensor 

These four requirements taken together 
appear to be just sufficient to restrict 
the overall form of the fundamental acoustic 
equations (i.e., the acoustic analogy 
equations) to trivial variations of the 
linearized hyperbolic conservation laws (6) 
to(8). 

Interpretation of the right side of these 
equations as acoustic sources implies that 
that the left sides account for all of the 
interaction of this sound with the base flow- 
including the propagation of the sound 
through that flow. The apparent linearity of 
the equations can be exploited to separate 
out these interaction/propagation effects 
from the unsteady source fluctuations that 
produce the sound by using the 4 th 
component of the adjoint vector Green’s 

functiong° V)a , which is related to 

the 4 th component of the ordinary vector 
Greens’ function g 4v (.*;,£ |j,t) by the 

reciprocity relation g a v4 

= g 4v and satisfies the 

inhomogeneous linear equation [9] 


4 



K\ t gU(yA x ’ t ) = 

(19) 

S M< S(x-y)s(t-v) 

with Delta-function source term where V 

is the adjoint of the operator L^ v that 

appears in equation 

Error! Reference source not found.. And 

setting aside solid surface interaction effects, 
for simplicity sake, we can use Green’s 
formula to show that the pressure like 
variable p' e is given by the tensor product 


Pe{x,t) = - { J y Vij (y,T\x,t)el J (y,T)dydT 

( 20 ) 


V —00 


of the source strength (13) with the 
propagator y^- which is related 

to the adjoint Green’s function 






fy,- 


dv„ 


( 21 ) 


( Y_1 )p JLga 44(^ T 


with v 4 = 0 ,and, therefore, accounts for all 
of the base flow interaction effects. 

Notice that p' e ( x, t j reduces to the ordinary 

pressure fluctuation p'(x, t) when the 
observation point x is in the far field. 


3 Specific Analogies 

The present results are, of course, 
general enough to allow for many possible 
analogies— each of which corresponds to a 
different choice of base flow— with the most 


rudimentary being obtained by setting the 
base flow stresses 

e v j and the velocity v. to zero. The 

hyperbolic conservation laws then imply 
that the remaining base flow dependent 
variables p, p e will be equal to constants, 

say to their values at infinity. This is 
essentially the original Lighthill [1] 
approach as modified by Lilley [10,1 1] to 
separate out the isentropic and non- 
isentropic components of the 
pressure/density term that appears in 
Lighthill’s source function. This turns out to 
be important because it is only the latter 
component that can be associated with 
viscous and heat conduction effects that are 
almost universally neglected at the high 
Reynolds numbers of interest in most 
Aerodynamic noise problems. The former 
can be an important noise source that has to 
be modeled in this type of acoustic analogy. 

But the main criticism of the Lighthill 
approach is that the residual velocities v' 

that appear in the residual equation source 
term have steady flow components that 
introduce steady and linear terms into the 
residual stress tensor. But the former cannot 
generate any sound and therefore do not 
belong in the source function, while the 
latter are at least in part associated with 
propagation effects and might, therefore, be 
best book kept on the left hand side of the 
equations. 

This criticism is at least partially 
overcome by the next most complicated 
form of the analogy which amounts to again 
setting the stress tensor e vj - to zero but now 

requiring that the base flow dependent 
variables v. , p, p e satisfy the resulting 

Euler’s equations. However, the interest is 
usually in the steady flow solutions to these 
equations because for many turbulent flows 
the turbulent stresses act only slowly over 
long distances and the Euler equation 
solutions are, therefore, expected to provide 
a good local approximation to the actual 
mean flow. This is especially true for the 
important class of nearly parallel flows, such 
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as jets wakes and shear layers. In which case 
the relevant Euler equation solution is the 
uni-directional transversely sheared flow 

v i =S ij U(x T ),p = p(x T ) 

, p e = constant 

where subscript T denotes transverse 
components, the velocity v. = S^U (jc 7 , ) is 
in a single direction but, along with the 
density /? = /?( x T ) , can vary in planes 

peipendicular to that direction. The 
hyperbolic conservation laws (or in this case 
the Euler equations) then imply that the 
pressure like variable p e will be constant. 

The linear operator on the left hand 
side of the residual equations now reduces to 
the usual Rayleigh operator of linear 
stability theory fame and, as is well known, 
the residual equations can be reduced to a 
single inhomogeneous equation for the 
pressure like variable p' e which, in the 

Aeroacoustics context, is usually referred to 
as a Lilley’s equation [10]. And while it can 
be argued that this equation has served as 
the basis of most of the theoretical jet noise 
research for the past two decades, it has a 
number of issues which will be discussed 
subsequently. But it is first necessary to 
finish cataloguing the various specific 
analogies. 

The true reduced order modeling 
comes into play when specific turbulence 
models are introduced for the base flow 
stress tensor. As noted in the introduction, 
these models are usually based on the 
filtered Navier Stokes equations, which are 
obtained by applying a linear filter to these 
equations. From our point of view a linear 
filter is nothing more than a linear 
transformation of the continuous functions 
into a simpler subspace that commutes with 
differentiation. The resulting filtered Navier- 
Stokes equations are special cases of the 
general base flow equations with the stress 
tensor now given explicitly by 


etj =-p 

v i v j -v i v j ) + 




n -q < 23> 

V -V 

1 




-(y-1) 

"Ol 

Tl 

1 

1 ~ ~ ~ _1 
+ ^ e n V j+ e j V i 



(24) 


and with p e related to p by 

P e =P + ^Y~p[ v2 -^) (25) 

where the overbars denote filtered variables, 
the tildes denote the Favre (ref.2) filtered 
variables 

• = (P*)/P ( 26 ) 

and the usually unimportant viscous terms 
have, for simplicity, been omitted. 

It is possible to choose a filter that only 
filters out the components of the motion that 
actually radiate to the far field [12]. So the 
base flow will be completely silent and all of 
the residual flow will radiate away as sound. 
In which case, it would not be completely 
unreasonable to interpret the apparent source 
terms on the right side of the residual 
equations as the “true sources of sound”. 

But no one has as yet succeeded in 
developing a base flow closure model for 
this type of filter. While equations (24) and 
(24) are very suggestive, they do not 
actually provide a means for constructing 
such a model — beyond, perhaps, providing 
some constraints that they have to satisfy. 


4 Actual Mean Flow Analogy 

There is much more that can be done with 
this. But I’m definitely not going to get into 
the vast subject of turbulence modelling 
here, except to point out that the most well 
known filter, namely the time average 

1 T 

( 27 > 
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which is the quantity that is usually 
measured in Aeroacoustic experiments, can 
be expressed as the convolution product 

OO 


also produces a completely silent base flow, 
but it also filters out much of the non- 
radiating (i.e., non-acoustic) components of 
the motion. The effective pressure is now 
given by 

Pe=P' + ^-( pv' 2 -p^) (28) 

and more importantly, the base flow stress 
tensor e . is now equal to the time average 

e' y j of the generalized residual 

velocity/enthalpy Reynolds stress. So the 
actual stress tensor (acoustic source 
strength) that appears in the residual 
equations (6) to (8) is just the difference 


p 2 (x,t)= J 

—00 V 

(. y;tf,r)dy dq dr 

(31) 

of a “propagator” y v/|j/ , which is related to 
the original propagator y . by the 
convolution 

00 

J Yv/(*|M +t + ^)y li i(x\y + n,t l )dt l 



(29) 


(32) 


between this Reynolds stress and its time 
average— which means that it has zero time 
average as would be expected for a true 
acoustic source. 

The base flow equations now turn out 
to be the usual Reynolds averaged Navier- 
Stokes equations for which the turbulence 
modelling is certainly highly developed [3]. 
This form of the analogy completely 
removes the mean flow propagation effects 
from the source term, which is expected to 
significantly improve the medium and lower 
frequency predictions, but a steady base 
flow analogy cannot provide the requisite 
structure for inputting the long range 
turbulent scattering effects, which tend to 
become more important at the higher 
frequencies 

The Greens’ function solution (20) can 
now be used to show that the far field 
pressure autocovariance 
j t 

P 2 {x,t o) = — J p'(x,t)p'(x,t + t 0 )dt, 

11 -T 


(30) 


with the two point time-delayed correlation 

1 T 

*y/h / Cp; *) = — f <j (. y > K (? + n. h + 0 dx 0 

(33) 

of the source strength e" tJ ( y,z) , which is 

directly related to the generalized velocity 
/enthalpy Reynolds stress autocovariance 
tensor 


j T 

Kjni ^ T ) = ^7 1 [p v >'j - p v >'j ]( y> T o ) 


-T 

r , f 


P V, -PV/ 

by the relatively simple linear transform 


(y + v\,z 0 +T)dT 0 

(34) 


*vyp/ ^v/p/ 

2 ( 8v /%p/ + ( \i/ ^v/kk 


J 5 y/ 5 p l R iikk- 

(35) 
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Since this autocovariance tensor is about as 
close as you can get to what is actually 
measured in turbulent flows, these results 
provide an essentially exact relation between 
the quantities that are typically measured in 
Aeroacoustic experiments — as would be 
expected from any good theory. 

The tensor R v ■ ^ ( y; t], r) is also the 

quantity that ultimately has to be modeled in 
the acoustic analogy approach. The hope is 
that most of the non-local “propagation 
effects” which would be very difficult to 
distinguish from the turbulent fluctuations— 
and, therefore, very difficult to model- have 
been removed from these stresses. 

But the result (3 1) to (33) is actually 
more general than this and applies to any 
steady base flow, including the uni- 
directional mean flow. This greatly 
simplifies the computation of the 
propagator y . z , but, of course, only applies 

to nearly parallel shear flows. And even in 
that, admittedly important, case has, as 
mentioned in the previous section, a number 
of serious issues. First of all the source term 
contains steady as well mean flow 
interaction terms because v' still has a mean 


(steady) flow component. A more serious 
difficulty is that the Rayleigh operator that 
now appears on the left hand side of the 
residual equations can support linear 
instability waves that grow without bound 
far downstream in the flow. And when the 
Greens’ function is required to satisfy 
causality, which seems like a reasonable 
thing to do, the corresponding homogeneous 
solutions to the Rayleigh equation contribute 
to the adjoint Greens’ function and therefore 


to the propagator y (y, j causing 

them to become infinite there as well. An 
even more serious difficulty at supersonic 
speeds is that the Rayleigh operator has a 
singularity at the so called critical layer 
which produces a much stronger non- 
integrable singularity in the propagator 



xj 'j when the observation point 


x is in the far field — leading to the 


ridiculous conclusion that the radiated sound 
is infinitely loud. 

But none of these difficulties would 
occur it the base flow were taken to be the 
actual mean flow field — even in the 
important case of a nearly parallel shear 
flow. Because, while the linear instabilities 
can initially grow in such flows, the slow 
divergence of the flow will eventually cause 
them to decay — keeping the Green’s 
function and, therefore, the propagator 
finite. And even more importantly at 
supersonic speeds, these flows do not posses 
critical layers and, therefore, have no critical 
layer singularities. But the propagator is 
much more expensive to compute in this 
case — in fact, prohibitively so for many 
applications. 

5Application to Slowly Diverging 
Mean Flows 

There is, however, a way to obtain the 
best of both worlds — which amounts to 
using a perturbation approach that takes 
advantage of the small spread rate, say £ , of 
the mean flow. The base flow velocities 
will then expand like [7,8] 


v, =U(Y,y T ) + sU (x) (Y,y T ) + ... (36) 


v r = sV (Y, y T ) + £ 2 V (1) (Y,y r ) + ... (37) 
where 

Y = s y x (38) 

denotes a slow streamwise (source) variable 
that varies on the streamwise length scale of 
the mean flow , and 


yr={y2,y 3 } ( 39 ) 

v r = {v 2 ,v 3 } (40) 

denote cross flow variables. Of course, the 
other base flow dependent variables will 
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have similar expansions. And it is assumed 
that all lengths have been normalized by 
some characteristic cross flow dimension 
and that all velocities have been normalized 
by some appropriate characteristic 
streamwise velocity with similar obvious 
normalizations for the density, pressure and 
temperature. The lowest order terms in this 
expansion correspond to the unidirectional 
transversely sheared flow with the slow 
streamwise variable entering only 
parametrically. 

It might be expected that the 4 th 
component of the adjoint Greens’ function 
would posses a corresponding expansion of 
the form 

gy M =gva+^gv.+- ( 41 ) 

But in so far as this expansion is concerned 
the lowest order base flow solution still acts 
like a uni-directional transversely sheared 
mean flow. So nothing much has been 
accomplished here because the lowest order 
term satisfies Rayleigh’s equation and, 
therefore, still posses a critical layer 
singularity. But embedding this solution in 
this more global context allows us to 
construct a new “inner solution” in the 
vicinity of the critical layer that brings in 
nonparallel mean flow effects to eliminate 
the singularity there. Standard singular 
perturbation techniques [13] can then be 
used to combine these two solutions into a 
single uniformly valid result that remains 
finite everywhere in the flow and is not that 
much more expensive to compute than the 
original parallel flow result. The resulting 
“propagator” still becomes large when the 
source point is at the critical layer but, 
unlike the parallel flow result, now remains 
finite there. 

This result can now be inserted into the 
formula (31) for the pressure autocovariance 
which can then be Fourier transformed to 
show, upon assuming only that the 
transverse length scale of the turbulence is 
short compared to the transverse length scale 
of the mean flow, that far field acoustic 


spectrum is given be the following purely 
algebraic expression[8] 

7 co(*W-> 


^ 2n^~ 
v X J 


sin e r Kj (* \y T J ) r L (* \y T J) 


0 * 


1C jXl 


y;—cosQ,—VS,(o(l-M c cos0) 

L c* c w J 


as x — » co (42) 

where S denotes the Eikonal of the 
corresponding geometric acoustics solution, 
the turbulence source correlations enter only 
through the spectral tensor 


1 00 

2n J l 


j{k&i+ k T< T ) -M 


(43) 


where the asterisk is being used to denote 
complex conjugates, r ■ is 

proportional the far field expansion of the 
Fourier transform of y • , x) is 

the moving frame correlation defined in 
terms of the fixed frame correlation (33) by 
[ 8 ] 

*%Ay>f> x ) s 

/ ' \ ( 44 ) 

^vjvd\y^ +iu c^) 

(where U c denotes the convection velocity 
of the turbulence) and I have introduced 
4 ( x | j) ,the acoustic spectrum at the 

observation point x due to a unit volume of 
turbulence at the source point y ,for 
transparency sake. Obviously, this quantity 
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has to be integrated over the entire noise 
producing region of the flow to calculate the 
actual acoustic spectrum 

= (45) 

v 

And finally 1 have neglected the contribution 
of the linear instability waves primarily 
because they do not seem to be all that 
important at the relatively low supersonic 
Mach numbers that we will be dealing with 
here. 

6 Application to the Jet Noise 
Problem 

The results of the previous section apply 
to any parallel shear flow but, for 
definiteness, the reminder of the paper will 
be restricted to the technologically important 
jet noise problem. The propagator 

( x | y T , Y ) ultimately depends only on 

the mean flow which, in today’s 
environment, would almost certainly be 
determined from a RANS computation. But 
the spectral tensor 

®y/u/ {y>ki’k T ,(o) depends on the 

turbulence statistics, which can not be 
obtained from any steady flow solution. It, 
therefore, has to be modeled, but the model 
can be paramatized and the parameters can 
be determined from some reduced order 
flow computation such as a RANS solution. 

Ideally, we would like to model the 
spectral tensor ® * ijkl itself, but the models 

must, at least at present, be based on 
experimental data and the experimentalists 
are unlikely to measure this quantity any 
time in the near future. A fall back position 
might be to develop models for the spectra 

^vt/U'Mr’Co) 

f +kT ’^R^{y,^)d^dz 

2n io r 

( 46 ) 


of the generalized Reynolds stress 
autocovariance tensor, which can be related 
to ® * t j kl by using the linear transform (35). 

But even this quantity has only been 
infrequently measured (see, for example, 
Harper-Bourne) [14] -and then only for very 
low Mach number flows. The only recourse 
is to model the Reynolds stress 
autocovariance tensor itself and to work 
backward through these formulas to 
calculate the spectrum ® * i - kl that actually 

appears in the acoustic equation(42). 

There are two main requirements for 
such models. The first is that they must 
reduce the large number of independent 
spectral components to a manageable level 
(there are 45 of these in all even when the 
enthalpy fluctuations are neglected, there are 
78 when they are not) and the second is that 
the models for the remaining components 
must be relatively inexpensive to compute. 

The first of these is usually met by 
introducing some sort of symmetry and/or 
statistical assumptions such as local isotropy 
and quasi-normality. But these assumptions 
do not seem to be all that viable at the high 
Mach number end of the range of interest. 
We, therefore, assume only that the 
turbulence anisotropy is Mach number 
independent and that the turbulence itself is 
axysymmetric [15] — which is consistent 
with the experimental observation that the 
turbulence statistics in the various cross flow 
directions are much more similar to one 
another than they are to those in the 
streamwise direction. 

The second requirement usually implies 
that the integrations that arise in the ® * i j kl 

calculation can be done in closed form, or 
nearly so, because of the relative expense in 
doing numerical integrations over and over 
again. But the models must also be flexible 
enough to provide an accurate representation 
of the turbulence structure because the 
strong streamwise retarded time variations 
cause the radiated sound field to be very 
sensitive to that structure. 

Appropriate models were introduced in 
reference [8]. They involve a large number 
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of adjustable parameters (actually infinitely 
many) which, as noted above, have to be 
determined from some reduced order flow 
computation — which in today’s world would 
again have to be a RANS computation. But 
this type of calculation can only provide 
enough information to determine a small 
number of these parameters. However, 
newer higher fidelity reduced order flow 
computations, such as URANS, VLES or 
even hybrid RANS/LES, are rapidly coming 
on line and these codes should be able to 
determine many more of these constants in 
the near future. 

But the results that will be presented in 
this paper are based on a pure RANS code, 
namely the NPARK WIND code, which is 
widely used in the Einited States [16]. This 
code was used to calculate the mean flow 
from cold (i.e., unheated) round jets with 
acoustic Mach numbers Mj = Uj / c x of 

0.50, 0.90, and 1.4, where U y is the jet exit 
velocity and c x is the speed of sound at 

infinity. The upstream conditions were 
specified in terms of the nozzle temperature 
and pressure ratios. The results were then 
used in equation (42) to calculate their far- 
field acoustic spectra on the arc xl Dj= 100 

and compared with jet noise measurements 
taken NASA Glenn SHJAR rig [17] with the 
same upstream conditions and with the 
atmospheric absoiption removed from the 
data in order to make the comparisons on a 
lossless basis. The results are shown in 
Fig.l. The lower curves are the acoustic 
spectra at 90° to the downstream axis and 
the upper curves are the spectra at 30° 

(which is close to the peak radiation 
direction). 

The overall agreement appears to be 
quite good but there is a tendency to under 
predict the higher frequency components of 
the 90° spectrum in the supersonic case. This 
is because the flow is not correctly expanded 
here and the present analysis does not 
account for the resulting shock associated 
noise (shown cross hatched in the figure). 

Computations were only carried out at the 
relatively low (<1.5) supersonic (acoustic) 


Mach numbers where non-linear 
propagation effects are believed to be 
unimportant. Fortunately, this also 
corresponds to the Mach number range of 
most technological interest. The critical 
layer only appears when the observation 
angle (as measured from the downstream 
axis) is fairly small (<45°) and gradually 
moves inboard from the nozzle lip line with 
increasing downstream distance until it 
reaches a point beyond the end of the 
potential core where it quickly moves onto 
the jet axis and suddenly disappears. But the 
propagator is still very close to being 
singular, and can consequently be relatively 
large, for a significant distance downstream 
of this point. It was, therefore, necessary to 
construct an additional inner solution for this 
region as well — especially since much of the 
small angle sound field will actually be 
generated in this relatively localized portion 
of the jet when the acoustic Mach number is 
close to unity. 

Finally, in order to put these results into 
perspective Fig. 2 shows a comparison with 
the best previous calculation using the Glenn 
JeNo code with an ad hoc correction to 
eliminate the critical layer singularity [18]. 
The results clearly demonstrate that there is 
enormous improvement in predictive 
capability with the present code. 
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Fig. 2) Comparison with JeNo Results. (Taken From [8]) 
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